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In the cited paper by the late J. Kurkijarvi and myself (KA) a molecular dynamics method is used to simulate the 
thermal agitation of the phase variable in a current carrying Josephson junction. Since the method is general and 
(^ may be useful in other contexts, it is briefly described with some details added. 

fSJ Consider a classical particle of mass M, position X, and momentum P performing a one-dimensional thermal 

» , Brownian motion in a potential U. The process may be described by the Langevin equations 

C^ ' dU 

X^P/M; P=-—-TjP + L{t), (1) 

O; ' where dots indicate derivatives with respect to the time t, and the symbols not yet defined are rj the average dissipation 

and L(t) the fluctuating force. It is known that to lead to thermal equilibrium at temperature T (in energy units) the 

,—1 fluctuating force must have the autocorrelation function {L{t)L{t')) = 2r]MT5{t — i'), where the brackets indicate an 

,^ . appropriate average. 

^ ' The basic idea of the method introduced in KA is to simulate the fluctuating force L{t) by random impulses 

^ describing elastic collisions with an ideal gas of light particles of mass tti, assumed to be always in equilibrium. 

G ' From the collision dynamics one calculates that random impulses 2pi drawn from a set distributed according to 

g[p) — (|p|/2toT) exp(— p^/2?TiT), — oo < p < oo, occurring at random time intervals ti from a set distributed 

according to f{t) — i/exp(— i/i), < i < oo, are required. With a given choice of tti <^ M, v must be chosen to 

!/3 ■ satisfy v = rjM/Am, this being the relation between the mean frequency of impact v and the damping constant rj. 

Although the algorithm outlined in the last paragraph follows so directly from the physically motivated model that 
it cannot be wrong, it is worth verifying that the autocorrelation function of L{t) so produced is indeed correct. 
In the model L{t) = Y.Zi '^PiK't - Ti), T, = Y,-^^ tj. It follows that 
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{L{t)L{t')) = (4p2) X {5{t - h) + Sit -h- h) + 5{t -h-t2- h) + 5{t ~ h - t2 - h - U) + . . .) x S{t - t'). (2) 



The fact that (piPj) = 0, for i j^ j, has been used to keep only diagonal terms in a double sum. In Eq. ([2]), the first 

term (the average squared impulse) is independent of its index and of time: it has been taken out of the time-average. 

Its value from the given distribution is 8mT. 

" ' ' The average of the sum of (5-functions is v. This can be seen both intuitively and formally. The intuitive argument 

J^ I is that, since the mean time interval between the (5-functions is j/~^, the integral of the average over a time N/v is N. 

I . The average can also be done exactly term by term. The first ternr is simply 

^^ ■ For the last term shown the average is 

(N ■ 
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AA = {5{t-ti-t2~h~U)) ^ I dt2 dt3 dtiv'^e-''^Q{t-t2-h~U) 

Jo Jo Jo 

dt2 I dh / dU v^e""' ^T\''^ ^^ ^^''*- ("^^ 
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■ - - Above, the first equality comes from using the (5-function to collapse the ti integral, the step function Q imposing the 
requirement that ii be positive. 

[The average A4 may also be evaluated using Fourier transforms and a contour integration: 



_ 2tt oj + iv 

= "'^'i^'"'''*'"=-^'^ ^ "" ^M'e-'^*.] (5) 

As illustrated, the successive terms in the sum of (5-functions, when averaged over the random time intervals, 
produce an infinite series for exp(i/i) which exactly compensates for and removes the decaying exponential in Eq. ([3])! 
The intuitive argument that the average of the sum of (5-functions is equal to v is thus formally verified. 



Inserting these evaluations into Eq. ([2]) one obtains 

{L{t)L{t')) = SmTv5{t - t') = 2r]MT5{t - t'). (6) 

Note that in the final result the mass m of the bath particles no longer appears, but that the temperature of the 
bath has been communicated to the Brownian particle. 
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